NACA TN 4347 


- 


CASE F1LF 


COPY 


NATIONAL ADVISORY COMMITTEE 
FOR AERONAUTICS 


TECHNICAL NOTE 4347 


A NONLINEAR THEORY FOR PREDICTING THE EFFECTS OF UNSTEADY 

LAMINAR, TURBULENT, OR TRANSITIONAL BOUNDARY LAYERS ON 

THE ATTENUATION OF SHOCK WAVES IN A SHOCK TUBE 

WITH EXPERIMENTAL COMPARISON 

By Robert L. Trimpi and Nathaniel B. Cohen 

Langley Aeronautical Laboratory 
Langley Field, Va. 



Washington 
September 1958 


naca re copy 


Loan ev 

4, rn _ 


PEASE®' 1 * J 


won 

HE 


ru 


Washington 


ll, 0. C. 


NATIONAL ADVISORY COMMITTEE FOR AERONAUTICS 


TECHNICAL NOTE 

A NONLINEAR THEORY FOR PREDICTING THE EFFECTS OF UNSTEADY 
LAMINAR, TURBULENT, OR TRANSITIONAL BOUNDARY LAYERS ON 
THE ATTENUATION OF SHOCK WAVES IN A SHOCK TUBE 
WITH EXPERIMENTAL COMPARISON 
By Robert L. Trimpi and Nathaniel B. Cohen 

SUMMARY 


The linearized attenuation theory of NACA Technical Note 3375 is 
modified in the following manner: (a) an unsteady compressible local 

skin-friction coefficient is employed rather than the equivalent steady- 
flow incompressible coefficient; (b) a nonlinear approach is used to 
permit application of the theory to large attenuations; and (c) transi- 
tion effects are considered. Curves are presented for predicting atten- 
uation for shock pressure ratios up to 20 and a range of shock-tube 
Reynolds numbers. Comparison of theory and experimental data for shock- 
wave strengths between 1.5 and 10 over a wide range of Reynolds numbers 
shows good agreement with the nonlinear theory evaluated for a transi- 
tion Reynolds number of 2.5 x 1C)6. 


INTRODUCTION 


The increasingly widespread use of the shock tube as an aerodynamic 
testing facility has led to the closer investigation of the flows present 
in such tubes. In particular, since the deviation of these flows from 
those predicted by perfect fluid theory is often of large magnitude, 
these deviations have been investigated fairly thoroughly. Several such 
studies, either of an experimental or theoretical nature, may be found 
in references 1 to 11. Investigations of the boundary layers in shock 
tubes have been made in some of the aforementioned references as well 
as in references 12 to 17 . This list of references does not cover the 
complete field of literature existing on these topics but is representa- 
tive of the various general treatments. 

Consideration of the entire flow field from the leading edge of 
the expansion wave to the shock wave is necessary to obtain an accurate 
picture of the waves traveling along the shock tube. These waves are 
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responsible for the deviations from perfect fluid flow in shock-wave 
strength (attenuation) with distance, in pressure ahd density at a given 
distance with time, and so forth, which have been noted by various 
investigators. The analysis of reference 1 was the first to treat this 
complete flow. Figure 1, a reproduction with minor changes of figure 1 
of reference 1, is the basic wave diagram of the unperturbed shock-tube 
flow showing the various flow regions to be considered with typical char- 
acteristics and particle paths. This linearized analysis (ref. l) was 
based on an averaged one-dimensional nonsteady flow in which wall-shear 
and heat-transfer effects generated pressure waves to perturb the' perfect 
fluid flow. This averaging process essentially implied thick boundary 
layers. The expansion wave was treated as a "negative shock" or zero- 
thickness wave. The resulting perturbation equations then hinged on the 
evaluation of the local skin-friction coefficient cp, which in refer- 
ence 1 was assumed to be given by an equivalent incompressible steady flow. 
Consequently, the application of the results of reference 1 was limited 
to shock pressure ratios in which this assumption for Cf was valid, 
although the analysis was still applicable for other pressure ratios 
when the proper choice of cp was employed. The assumption of incom- 
pressibility should apparently eliminate the strong shock pressure ratios 
from the range of validity. 

Solutions to the laminar boundary-layer equations employing a linear 
viscosity-temperature relation (refs. 2, 3 > , and 14) show that the non- 
steady character of the flow is such that the equivalent laminar steady- 
flow assumption is in error, irrespective of compressibility, for most 
conditions except that existing in the cold-gas region a for strong shock 
waves. On the other hand, the turbulent boundary layer is not nearly so 
sensitive to the unsteady character of the flow. Reference 15, which 
assumed a one-seventh-power velocity profile similar to that of refer- 
ence 1, reported that even for infinite shock pressure ratios the effect 
of unsteadiness would produce only a maximum variation in turbulent skin 
friction of 5 percent in the cold gas and of 22 percent in the hot gas. 

The o nl y other attenuation analysis to date that considers the entire 
flow field is that of reference 2 . This analysis is similar to that of 
reference 1 in that it is a small-perturbation approach using traveling 
waves and a negative shock, the major difference being that the pressure- 
wave generations arise because the boundary-layer-displacement thickness 
changes with time. (The boundary-layer-displacement thicknesses of 
ref. 15 are used.) Flows with thin boundary layers having a linear 
viscosity -temperature variation are required for this treatment to apply. 
The attenuations predicted by references 1 and 2 for turbulent boundary 
layers agreed within 10 percent for shock-pressure ratios up to 6 in 
spite of the marked differences assumed in the mechanism for handling 
the wall effects. The perturbations in the flow behind the shock show 
a larger difference between the two approaches. 
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The deviations from ideal theory discussed arise for the most part 
from vail effects, that is, the perturbations in the shock-tube flow 
caused by vail shear and heat transfer. Much recent vork has been done 
using the shock tube as a testing medium to provide very high-temperature 
flovs of short duration. (See, for example, refs. 9 en<l l8«) In these 
cases, deviations from ideal fluid flov vill also arise because the air 
at high temperatures does not behave as an ideal fluid. It vould be 
difficult to separate the real-gas effects from the vail effects; there- 
fore, the present analysis, like those of references 1 and 2, is con- 
cerned only vith the effects of vail boundary layer upon the inviscid 
outer flov, the fluid being considered as an ideal gas. 

The turbulent theory of reference 1 has been compared vith experi- 
mental data for attenuation in references 1, 7> and 10 and good agree- 
ment has been found in general. Predicted pressure perturbations in the 
hot gas by the method of reference 1 agreed veil vith the experimental 
results reported in the same paper. Fair to good agreement betveen 
theory and experiment is reported in references 7 and 10 for the hot- 
gas average density variation vith time in the flov behind the shock 
vave; poor agreement is reported for the cold-gas flov vhere the finite 
expansion fan has been treated as a negative shock. 

Since the deviations from the inviscid fluid flov often become 
large in cases of aerodynamic shock-tube testing, the linear, or small- 
perturbation, theories of references 1 and 2 are no longer applicable 
and recourse must be made to some sort of nonlinear approach. 

In order to obtain an exact theory for predicting the perturbations 
in a shock-tube flov, a rigorous treatment vould be required first to 
the solution of the boundary-layer flovs. The boundary-layer equations 
vould have to be solved not only in region 3 but also inside and after 
the expansion fan vhich is considered to be of finite extent. For 
laminar flovs the main difficulty vould probably be the correct handling 
of the viscosity variation across the boundary layer. For turbulent flov 
a rigorous treatment appears to be impossible vithout a tremendous 
increase in knovledge of the mechanics of turbulence. Once the boundary- 
layer solutions vere determined, the vertical velocity at the edge of the 
boundary layer could, if the boundary layers vere thin, be used in the 
manner of reference 2 to determine the local pressure vaves generated. 

The second major difficulty in obtaining a rigorous perturbation 
solution vould arise from the treatment of the entropy discontinuity. 

The theoretical contact surface increases in extent vith distance pro- 
gressed dovn the shock tube due to mixing and diffusion (the former is 
the major influence). This process not only generates pressure vaves 
but also alters the reflected and transmitted vave strengths of the pres- 
sure vaves generated by the boundary layers. 
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If a rigorous solution such as that just described was available, 
then it would without question be the one to be employed. The theory of 
reference 1 and the theory of the present report assume that the wall 
effects can be averaged across the flow. This assumption introduces errors 
because no such physical mechanism exists for the instantaneous trans- 
mission of these effects across the flow. In the absence of the rigorous 
solution there is no evidence to indicate that the errors introduced by 
a shock-tube perturbation theory based on an averaging process are of a 
larger magnitude than those introduced by the neglect of the aforementioned 
considerations required for a rigorous solution. In addition, there is 
the possibility that an averaging process might be more applicable as the 
boundary layer fills a greater part of the shock-tube cross-sectional 
area. Consequently, the extension of the method of reference 1 in the 
present analysis is justified. 

In the present paper the analysis of reference 1 is first modified 
to eliminate the restrictions imposed by the incompressible equivalent 
steady-flow assumption for local skin friction; and then a nonlinear 
theory is derived which permits application of the analysis to large 
attenuations. It will be assumed that the reader is familiar with the 
basic theory and assumptions of reference 1 so that repetition in this 
paper may be avoided. This modified theory will be compared with experi- 
mental data covering a wide range of flow variables. The theoretical 
and experimental studies reported herein were conducted at the Gas 
Dynamics Branch of the Langley Laboratory during 1955 and 1956. 


SYMBOLS 


a 



Cn 

c f 

c v 


c 


P 


velocity of sound 


constant defined by equation (8) 

local skin-friction coefficient, 2-r w /pU^ 

coefficient of specific heat at constant 
volume 

coefficient of specific heat at constant 
pressure 
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D 


hydraulic diameter. 


4 X Area 
Perimeter 


E = 1 + 


. XI 


ify - u\e 

P Xn> §Xm> §Xn 


Fn = 2 


*— r n+ b 
n + 1 e^p S 3 

nT~3 b) ° n 


G m or Gn 

Sn = K S n 
K 


Ln or T,n or T 


l 

l* 


functions defined by equation (39b) 


constant defined by equation ( 36 ) 


ratio of contributions of P waves to total 
waves generated in region 0, 

(M s - Mg /0p,p \ 

1 + M p - (M s y Pp J 

linear attenuation with first subscript 
describing boundary layer appropriate to 
region a and second subscript to region 0; 

that is, t = Pys " P P° for region a 

^OO 

turbulent with n = 7 and region 0 as 
transitional 

fixed distance along shock-tube axis 

distance the shock moves from a given point 
until the effects of transition in the flow 
generated at that point first influences 
the shock 


M = U/a 
Ma = U/aa 
Mp = U/ap 
Ms = bg/ &00 
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m or n 

reciprocal of velocity exponent in boundary 

layer> 5 = (l) m 

N 

arbitrarily denotes subdivisions I, II, III, 
etc . of hot-gas region 

NLn or T,n or T 

nonlinear attenuation with first subscript 
describing boundary layer appropriate to 
region a and second subscript to region P 

N Pr 

Prandtl number 

P,Q 

2c v 

characteristic parameter, a ± U 

R 

A A 

p,Q 

effective characteristic wave parameter 
defined by equations (43) 

p 

static pressure 

R 

gas constant: Reynolds number 

V 

R* 

Reynolds number of transition, 

S 

entropy 

T 

temperature 

T v 

wall temperature 

v 

adiabatic wall temperature 

t 

time 

U 

free-stream velocity 

u s 

shock velocity 

U W 

velocity of wave which generates flow 

u 

X 

velocity in boundary layer 

distance along shock tube from diaphragm 
station 


x 
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y 

Z N 

r 

7 

a 1 


distance from surface 

function defined by equation (B12) 

function defined in equations ( 31 ) 

ratio of specific heats , Cp/c v ; assumed 
as 1.40 for computations 

length of segment into which shock tube is 
divided for nonlinear treatment 


A( ) 
AX 


A( ) 

A| 

6 

S* 


S( ) 

5t 



characteristic derivative in boundary layer. 

) , 1 &* a( ) 

ax u e at 


characteristic derivative in boundary layer, 

A. + i 1 d( ) 
ai u e e at 

boundary- layer thickness; also indicates 
differential quantity 

boundary- layer displacement thickness, 



characteristic derivative in potential flow, 

<0 + (u ± a)^-A 
St Sx 


A 

contribution to attenuation due to P waves 
in region (3 


1 



nr 


recovery factor, assumed equal to O .85 for 
laminar flow and O. 9 O for turbulent flow 
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0 


9 0 


boundary- layer momentum thickness, 



boundary-layer momentum thickness at | 


k = — for £* > ^ and if £* ^ its value is 1 


d 

v 

I 

^d 

^o 

I* 

P 


coefficient of viscosity 

coefficient of kinematic viscosity 

distance flow has progressed along surface 

distance flow has progressed along surface 
at entropy discontinuity 

distance flow has progressed along surface 
when 0 = 0 O 

distance flow has progressed along surface 
at transition 

density 


a = — 
a p 

T W 

0P,oi>0P,|3>0Q,P’0s,a>0s,|i 
tyn or 


vail shearing stress 

influence coefficients, defined in appendix A 

influence coefficients defined by- 
equation (44) 

compressibility correction 


a) exponent in viscosity-temperature law, 

H « 

subdivisions of hot-gas region 0 for non- 
linear treatment (see fig. 8) 


I, II, III, etc. 
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Subscripts: 

Letter subscripts not included in the symbols defined above refer, 
in general, to values at points or within regions shown in figure 1. 
Numbered subscripts refer to points in figures 8 and 9* Exceptions to 
be noted, however, are as follows: 

m ^ n refers to velocity profile parameter m,n 

o perfect -fluid value 

-t at time t 

y g evaluated immediately behind shock, that is, 

point v located at x = U s t 


x at distance x 

1 denotes attenuation with transition 

X arbitrary condition in shock-tube free stream 

g-j-^ denotes NACA standard atmospheric conditions 

A prime on a symbol indicates a quantity evaluated at reference 
temperature . 


THEORY 


Derivation of Expressions For Local 
Skin-Friction Coefficient 


The skin-friction coefficient for the flow behind wave-induced flows 
be found by an integral method. An incompressible skin-friction 
coefficient will first be determined and then a simple compressibility 
correction will be applied. 

The integrated equation of motion for the incompressible boundary 
layer with zero pressure gradient is (see ref. l4): 

(M'l + itel) . In- = i c f 

\8x/ t u\8t / x pU 2 2 


( 1 ) 
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The form parameter 6*/0 is assumed to be constant; this has been 
shown to be true for the unsteady wave-induced laminar flow (see ref. l4) 
but has not been completely established for turbulent flows. Equation (l) 
then becomes 



( 2 ) 


Since the resulting expression for cp will ultimately be used in 


the attenuation formulas wherein the integral 



is desired. 


the variable | is introduced. The variable £ is defined as the dis- 
tance a particle in the free stream has moved to reach the point (x,t) 
since acceleration by the passing wave which originated at x = 0 at 
t = 0 and which travels with velocity Ify. Thus, 


l(x,t) 


— (Uyt - x) 

UW - U 


( 3 ) 


In the case of flow in a shock tube, the value of Uw is U s for 
the flows induced by shock waves. If the assumption of reference 1 is 
followed and the expansion wave replaced by a wave of zero thickness, 
moving with the speed of the leading edge of the original expansion 
wave, then tfyj = -a e . 


The differential equation (2) is transformed from the x,t coor- 
dinate system to the £,t system by using the following derivatives: 



(Sr') 

W Max 


U 

Uw - u 





. (ah (M'l t (V) . + (— H 

x Wt\dt/ X Wg % - u\^/ t 4 

> 

Thus equation (2) becomes 

fa?'; 1 1 5 ! i/m) , n c . 

yd£/ t u e E^dt/^ 2E 1 



(*> 


( 5 ) 
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where 


E = 1 


+ ^ fb* 

+ % - u \e 


(6) 


Now equations (2) and (5) are differential equations capable of 
solution by application of the method of characteristics. The slopes 
of the characteristics are 


At _ 1 6^ 
Ax U 9 


for equation (2) and 

At .1^1 
At U e E 

for equation (5)- 

Thus, if the symbol Af/A| is used to denote the derivative along 
a characteristic of slope At/A£, equation (5) ma y written as 


A0 = _1_ 
A| 2E 


(7) 


For steady-flow boundary layers with zero pressure gradient, it has 
been established that Cf = cp(6,U,v). If this relation is assumed to 
hold for unsteady flows in the same form as for steady flows, then 

< 8 > 

For turbulent flows n is the reciprocal of the exponent in the 

fractional power expression ^ used to describe the turbulent 

boundary-layer velocity profiles. For laminar flows the value of n is 
one. The C n terms are arbitrary constants to be evaluated later and 
may be completely different for the steady and unsteady flows. 
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Combining equations (8) and (7) and integrating yields: 

_ JL j|>_ 

A| 2 E \v 0 / 


( 9 ) 


e^A e.Sifvep r* 

253 V u8 / 


(10) 


0 = 


n +3 

n+1 | n + 3 ^n/v 0 
+ n+1 2Em 5 


.n+1 


(l - to) 


n+1 

n +3 


(11) 


Substituting equation (ll) into equation (8) to obtain c f as a func- 
tion of ^ yields: 


c-(S) = 2 


— t n+1 

2_L1 E «f + c n n+3 

n +3 8 n 


n +3 


n +3 


n + 1 a n+1 2E /5\ n+1 ^u\ n+1 u, 

n~+l 9o de v + ^ " *■ 


2 

n +3 


c f (|) = FnE^ 


2 p+ 3 

n + 1 2E/5\ n+1 ^U0 o ^ n+1 u A x 

(— J + v( 5 ‘ 5 o) 


2 

n +3 


(12a) 


(12b) 


For the special case of 9 0 = 0 at | Q which corresponds to flow 
initiation at | Q there results 


c f (|) = F n E n+ 3 - 4 


2 

n +3 


(12c) 


The values of the various F n terms, which are directly related to 
the hitherto arbitrary C n values, will now be determined to match known 
solutions for certain limiting cases. If the ratio tfy/u becomes 
infinitely large, the solution must be the same as that for an infinite 
flat plate in contact with a fluid impulsively started from rest at time 
t = 0. 
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Under these conditions 

I - 60 = Ut 



The solution, known as the Rayleigh solution, to such an impulsive laminar 
flow over a plate is 


c p ( t ) 



(iM 



Rayleigh 


2.469 

/ 


Consequently, in order to match the Rayleigh solution for the laminar 
casq (n = l) 


F 


1 



1.128 


or 


F x = 0.718 


(15) 


On the other hand, if lfy/U becomes infinitely small, the solution 
must be the same as that over a semi-infinite flat plate in steady flow 
(that is, the so-called Blasius problem). For these conditions 

E = 1.0 


c f (6) = F : 


I - Sr 


2 

n+ 3 


n\ 


(16) 
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whereas the corresponding form of the Blasius solution for laminar flow 
is 




c-p = 0.664 

(-) 

\ ® /Blasius 


jK‘ - to) 

= 2.605 


1 

2 


(17) 


Thus, if the Blasius solution were to he matched, it would be nec- 
essary for 


F x = 0.664 

Two possible solutions are then available for the laminar 
incompressible-flow case depending upon which limiting value is matched: 

/% \ 

Rayleigh limit l — » 00 J : 


( c f ) n ,l(l) = 0.718(1 ♦ 1.169 ^ j ) H({ - to 


\ 2 



V 

U/. 

1 

1 

V 

[v( ? 


1 

2 


(18a) 


Uu 

Blasius limit > 0 | : 

V u 


( c f)n=l^) = 0.66411 + 1.605 


% 


1 
^2r 


% - U/ \v 


S(* - ••) 


In figure 2 the values for c^ 




1 

2 


(18b) 


as determined by four 


different means are plotted against pressure ratio across a shock in an 
air-air shock tube. The upper branch of the curves applies to region 3 
behind the shock, and the lower branch applies to region a behind the zero 
thickness expansion wave associated with the shock of strength Pp^p^. 
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In addition to the values determined from equation (l8) , there are shown 
values which would be obtained, if the fluid were assumed to be incom- 
pressible, by the integral method of reference 14 and the numerical solu- 
tion to the Prandtl boundary-layer equations. (The results of ref. 15 
are applicable to this numerical solution.) The agreement between both 
the curves of equation (l8) and the referenced curves is very good. How- 
ever, since the curve based on the Rayleigh limit gives a better approxi- 
mation in region (3, which will be shown to dominate the attenuation equa- 
tions, the Rayleigh values of Fp = 0.718 and 5*/0 = 2.469 shall be 
used for the remainder of this paper. Reference 17 also employed a nor- 
malized Rayleigh velocity distribution in the treatment of flow induced 
by shock waves. 


For the turbulent case two analogous limiting processes are not 
available in order to determine the values of F n . The turbulent boundary- 
layer theory is semiempirical and relies on experiments to supply con- 
stants for the resulting equations. Since no "Rayleigh- type" experiments 

% 

have been performed, there is no limiting process — ->co to apply for 
the turbulent case. There is, however, the semi-infinite flat-plate 


solution corresponding to the limit »0. 

U 


This solution (see refs. 1 


and 19) assumes a velocity profile — 
friction coefficient expressed as 



c f U) 


0.0581 



V 



and results 

1 
' 5 


in a skin- 


(19) 


The combination of equations (l6) and (19) results in a value of 
Fy = 0.0581. 

\ l/n 

it may be shown that 


9 = n 

5 (n + l) (n + 2 ) 


For profiles of the family 


U U 


S* m 1 
8 n + 1 

5^ = n + 2 
0 n 


( 20 ) 
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Therefore, the expression for the one- seventh-power turbulent veloc- 
ity profile skin-friction coefficient in the nonsteady incompressible 
flow becomes 


c f U) 


0.0581 1 


+ 


2 % \ 
7 % - uy 


1 

5 




(21) 


Since there was only a minor difference between the unsteady- flow 
values of cp for laminar flow when based on the limiting cases of 

Uy Uy 

— -» 0 and — — > 00, it is expected that the turbulent -flow agreement 
U U 

Uy 

would be just as favorable if results were available for — -too. Con- 

U 

sequently, equation (2l) is assumed to be a fairly close approximation 
to the correct answer. Equation (65) of reference 15 is very similar 
to equation (2l) but was derived in a different way. 


The skin-friction-coefficient relation of equation (19) corresponding 
to the one- seventh-power profile law is no longer valid at arbitrary large 
values of Reynolds number in incompressible steady flow. Instead a loga- 
rithmic law is often used. (See ref. 19*) However, since a power profile 
is easily handled by these methods, the skin friction on a semi-infinite 
plate for these large Reynolds numbers is found to be closely approxi- 
mated by the relation 


c f 


= 0.0186 


4 - So) 


1 

8 


( 22 ) 


which is compatible with the relation g - (y/s) 1 / 13 . If for consistency 

it is further assumed that 6*/ 9 has the value 15/13 (the value for 
n = 13) , then the unsteady skin-friction coefficient would be given by 

1 _i 

8 


c f (|) = 0 . 0186 ( 




ICO 

:?l 

1 

C\J 

4 - 0 

13 Uy - ui 

V 


(23) 


No claim is advanced that a l/l3-power profile actually exists at these 
higher Reynolds numbers; it is only necessary for equation (22) to be 
valid for steady flow and the value of 5*/9 to be 15/13* 
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It was shown in reference l4 that the skin friction in region a is 
correctly given by the method of characteristics using the boundary con- 
dition 0=0 for |=0 at x = -a € t (that is, on the zero-thickness 
expansion wave) for that part of a where x < Ut9/8* and the boundary 
condition 9=0 for 5=0 at x = U s t (that is, on the shock wave) 
for that part of a where UtQ/5* ^ x ^ Ut. Sketch (a) shows a boundary- 
layer momentum characteristic. In order to determine the unsteady fric- 
tion coefficient at point c in the region in question, Ut9/6* % x % Ut, 
the boundary condition of 9=0 on the shock wave (point a) is correct. 
In the analysis which follows, the boundary condition 0=0 on the 
expansion wave (point b) is used instead. Friction coefficients are 
shown qualitatively in sketch (b). The solid curve represents the case 
with the boundary condition on the shock whereas the dashed curve repre- 
sents the case with the boundary condition on the expansion wave. In 
the region in question it is seen that the differences are not serious. 
Because of the relatively sma l l contribution of region a to attenuation 
as compared with the contribution of region 0 and also in the interest 
of simplicity, this error is neglected. 


Boundary-layer characteristic. 




t 


Sketch (a) Sketch (b) 

In order to handle the transitional flows which occur behind the 
waves in a shock tube, some approximation for cp in the transition 

region must be employed. Any of a number of assumptions is possible in 
this region. However, in view of the many assumptions already present 
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in the theory, the least complicating supposition will he employed in 
this case; namely, an instantaneous transition is assumed and the value 
of c-f in the turbulent region will be taken as the value which would 
be present had turbulent flow existed since the initiation of flow. In 
other words, at the transition point the local c^ changes discontinuously 

from the laminar to the turbulent value and the value £ 0 = 0 then applies 
in both laminar and turbulent regions. This assumption was used in the 
logarithmic transitional curve for steady flow of reference 19* Figure 3 
compares the integrated skin-friction coefficients of reference 19 with 
the curves obtained by the various power laws and the foregoing transi- 
tional assumption. The agreement appears to be very good. 


A simple compressibility correction will be based on the intermediate 
temperature or T' semiempirical method. This correction, expounded in 
reference 20 for laminar boundary layers and in reference 21 for turbulent 
boundary layers, assumes that the incompressible skin-friction relations 
apply to compressible flow if the properties of the compressible flow are 
taken as some intermediate value between the wall and free-stream values. 
Thus, if the relation 


c f ee 




2 

n+3 


(24) 


applies to a steady incompressible flow, then 


c 


f 



(25) 


will apply to a steady compressible flow for a certain choice of the 
primed state. The following values of the intermediate temperature T 1 
are given (see refs. 20 and 21 ): 

For laminar flow: 


— = 1 + 0.032M 2 + 0.58( — - l) 


( 26 a) 
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For turbulent flow: 


TV 

T 


1 + 0.035M 2 + 0.45 



( 26 b) 


It may be shown that for an arbitrary reference state, which may be 
either wall or free stream, equation ( 25 ) may be put in the form 


C-P = 


' w 


? u2 



2 

n+3 


n +1 2 


' T \ n+ 5/(j.'\ n+ 3 

n \~) 



(27a) 


(27b) 


If a temperature-viscosity relation of the form 



is assumed to apply, then the steady-flow compressibility correction % 
becomes 


n+l- 2 co 



( 28 ) 


It is assumed that this T' method is also applicable as a compres- 
sibility correction for the unsteady-flow skin-friction coefficients. 
Compressibility corrections in hot and cold gas ^ are plotted against 
shock pressure ratio in figure 4 for an air-air shock tube. The value 
co = 0.8 has been used to compute these curves. 

The results of this section may be summarized by the following 
expressions for skin-friction coefficient: 
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where 


; f in=l _ 
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( c f)n=7 
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k 5; 


fU 
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I j / 


(29) 


= 0.0186 1 + 


i c fjn=13 


2 % \ 8 /u 8 
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n+1 2 

M n+ 3 /ni \ n+3 


Evaluation of the Linearized Attenuation Expressions 

The basic linearized attenuation expressions derived in reference 1 
are summarized in appendix A. The following expressions result from the 
first identity and subsequent substitution and manipulation of the per- 
tinent equations from the appendix: 


^vs ~ P30 ^vs - Pf Pf - Pc ^c ~ 25 Pb “ 2p 0 

= -r t + 


(30) 


Pys ~ P 30 _ 1 

ae "D J 0 


a po 

-0p o — s s F P g. J Fao + 2Ull ~ P ££ 0 S B 

oM s - M p - 1 r{ *>* 0 M S - M 3 + 1 Pp 0 + Q 3o 7 /& ^ i2 a £ ys ’P 

iaao/ 


c f,pU) d| + 


1 r s d 
D Jo 


a po 


Mg + 1 p Po + Q(3 o _ 7-1 a oo P fio ~ P oo 
^ P » a M<x + t| + 1 Poo + Qpo 7 /an \2 a e 


0s, c 


- 1 


c f , a U) dft 


(31a) 
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■vs 


- P pO 


a 


e 


1 

D 



r 3 c f,P 





r^d 

Jo r “ C f,a « 


(3») 


Since 


Pvs “ Ppo y a e p vs “ Ppo 


P Pc 


2 a po a e 


( 52 ) 


(see ref. l) the following relation results 


p vs ~ P pp fy a e 1 
p po \ 2 a po D y 


-1 


r K ^p r 5 d 

J,m / c f,P,m + r P,n / * c f,p,n + 

'■'0 K^p 


= r, 


o K la P 5 d 

P a,m J ^ c f,a,m d £ + p a,n # c f,a,n d ^ 


^d 


(33) 

in which the subscript m applies in the laminar range 0 6 = p 

and the subscript n applies in the turbulent range £ a p < la^p < Id. - 
The value of k is defined as t&jz* for I* ^ and as 1 for ^ 1^. 

The total linearized attenuation is thus made up of the sum of the 
effects of regions a and p. This relation may be expressed mathematically 
as 


Pvs Ppo 
p Po 



( 34 ) 
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where 



Z ag 
2 a Po 



c f,a,m 


+ r 


a,n 



c f,a,n 



(35a) 




c f,(3,m 


d£ + r 


3,n 




or 


(35b) 


Now, with expressed in the form 


c f,nU) = F nEn n+5 ^ £ 


2 

n+3 y 


T_ 

T' 


n+l-2o> 

\ n+3 


(36a) 


c f,n(l) = G n fi n^ £ 


n+3 


integration of equations (35) results in the following equations, 
subscript X designates either region a or region 3*) 


( 36 b) 

(The 
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Rearrangement of terms of equation (37) produces 


( Pvs - Ppo\ 

2 ” 

\(*J\ n+3 

V p °° * X,m,n 

Aw J 


1-1 


7 a € ^3o 
^ a 3o Poo 


m+1 m+1 


2(m-n) 




n+1 


r G n n + 3, u vA w-3/6 d \“-5 


n+1 

-sr 


(58) 


p Bo 


X,m, n 


n+1 

re 


2 

n+3 


1-1 




m+1 

*\dh-3 


2(m-n) 

a D\1 


m 


(iW-3)(n+5) 


SXrAn 


1 - 


n+1 

nS-3 


(39a) 


where 


U X v e 


SXm ^Xml o 


2 

m+3 


a ~ V X; 


= z Pgo m_+ 3 r G (M 
2a poPeo m+1 ^ *^2 / 


bH-1 

m+3 


(39*) 


Note that f^ is a function only of the shock pressure ratio Pp 0 / P oo 

and the value of m for any region X. The term g^ requires in addition 

a temperature-viscosity relation which may be either an exponential type or 
some other form, such as the Sutherland equation. 

For the case of no transition when the flow is either completely 
laminar or completely turbulent, equation (39 a ) reduces to 


•vs 


p Po . 


CH-1 _ 2 

' A" 5 /y>\" mt3 
J)/ Vco / 


s Xm%m 


2k 
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or 



n +1 

i'A n+5 /V' 


g 

n+3 

SXn^Xn 


(40) 


The attenuation in this case for a given initial value of is 

obviously dependent only on the two parameters expressed as hydraulic 
diameters of shock-wave travel z/D and shock- tube Reynolds number 

aeoD/Voo- 


For the case with transition three parameters are required to 
describe the phenomenon at a given initial value of Pp 0 ^P ro ; the third 

parameter is the transition Reynolds number R* and enters into the 
analysis in the following manner: 


£ * = u£ v _™_ 2 v_ _z_ 
Id v a^D Z U v ro Id 


(41) 


When equation (4l) is substituted into equation (39a) , the linear 
transitional attenuation relation becomes: 


( Pvs ~ p (3o\ 

\ Po ° 'X,m,n 


n+1 


n-1 




nri-l 


nH-1 


n+1 




m+1 




n+1 


n+1 


U v a 




D v* 


n+1 

n+3 


n+1 


1 . (hi JL _L\ n+5 ( R *)rH-3/i f«P\ ' ( K )n+3 


(42) 


a 7 

where the term — — — — is a function of p 6 lp M . Thus, at a given 
U Voo i d P 0 / 66 


value of Pp 0 ^P OT ^he linearized attenuation is a function only of Z/D, 
SooDj 


Voo, and R* . 


The attenuation functions g^, the compressibility corrections 

and their products are presented in table I for shock pressure 

ratios from 1.0 to 20.0 and for m equal to 1, 7, and 13 . It may be 
seen by inspection of equation (40) that, for given shock tube Z/D, 
Reynolds number a. o pjv oa , and no transition, the attenuation contribution 
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of the region X is directly proportional to the product (ftg)^. In order 

to demonstrate graphically the behavior of these compressible attenuation 
functions, they have been plotted in figures 5(a) and 5(h). The magni- 
tude of the contribution of region 3 to attenuation increases monatomically 
with increasing pressure ratio, is always negative in sign, and thus tends 
to increase attenuation. On the other hand, for low pressure ratios the 
contribution of region a tends to increase attenuationj above a shock 
pressure ratio of about 5-9> the trend is reversed and region a contri- 
butes compression waves that tend to decrease attenuation. This reversal 
is discussed more fully in reference 1. 

When no transition is considered, the attenuation function for the 
entire flow field is found by adding the contributions of regions a and p. 
For the cases where the profile exponent m is the same in both regions, 
the total compressible attenuation function has been computed and is 
shown in table I and in figure 6 for values of m of 1, 7> and 13 . This 
function is given as 


Pvb “ p Po 



(^)a,m + (%)p,m 


The results of using the methods of references 1 and 2 are also shown 
in figures 6(a) and 6(b) for values of m of 1 and 7, respectively. 


For the laminar case the curve in figure 6(a) representing the 
method of reference 1 falls far below that of the present report, pri- 
marily because of the importance of the neglected unsteadiness effects 
as discussed in the introduction. The results of reference 2 are also 
below that of the present report (approximately 25 percent for shock 
pressure ratios from 4 to 10 ) and show better agreement at higher pres- 
sure ratios (only 10 percent below for a shock pressure ratio of 20). 

Agreement between the methods of references 1 and 2 and the present 
method is better for the turbulent case (m = 7; fig. 6(b)). The neglect 
of unsteadiness has a smaller effect upon the results of reference 1, 
although the effects of compressibility still give significant devia- 
tions at shock pressure ratios near 10. The present results and the 
results of reference 2 are in agreement within less than 15 percent for 
shock pressure ratios up to 10 and then diverge to a 20-percent varia- 
tion at a shock pressure ratio of 20. 


Figures 5 and 6 show that the cold gas contributes only a small 
part of the total attenuation for Pp 0 ^P OT <20. At Ppo^P^ = 20 the 

relative a contribution is larger than that for p^/p^ < 20; however, 
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it is only about 4 percent, 15 percent, and 25 percent of the total for 
m = 1, 7, and 13, respectively. 


Evaluation of Nonlinear Attenuation Expressions 

The expressions derived previously are based on a linearized or 
small perturbation analysis. However, for many conditions encountered 
in the shock tube, the attenuation is no longer small. In order to maxi- 
mize the tota.1 available experimental testing time, most experimental 
work is done at values of Z nearly equal to the total length of the 
low-pressure side of the shock tube. At these large values of Z the 
shock strength often has decayed markedly from its value for small Z. 
Consequently, relations for the attenuation under these conditions would 
be very desirable. 

An approximate method to obtain the attenuation for the cases where 
the small perturbation analysis is invalid will be described. First, 
consider parameters P and Q which are related to P and Q by 


P = P 


Q = Q 


a S 
7 R 

a S 
7 R 


(43) 


A A 

(The parameters P and Q of this report are identical to the param- 
eters P' and Q' of ref. 1.) When equation (43) is substituted into 
equation (60) of reference 1, the following equations result: 



1 _a_ 1 6 £ + J;_ 6 U 
7 a £ p 6 t “ a £ 6 t 


(44a) 
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A 8P' 



2 _a 
D a_ 


Uc f \ £7 - 1)M + i]m + Np r 5 Tv ^ T ^ V > 


(Wrt>) 



= i Uc 
D 



(4kc) 


/X A 

As discussed in reference 1, the changes in P and Q are evident 
in changes by wave motion of velocity and pressure but do not indicate 
the various changes in entropy. For example, if equation (44a) is solved 
for 6p/p, 


bp _ 7 f^/SP + 
p 2 a la e 



(^ 5 ) 


The value of P is associated with waves moving with the flow at 
a velocity of u + a whereas Q is associated with waves of the opposite 
family moving at a velocity u - a. Now, for the linear attenuation 
theory it is assumed that reflections at the shock wave and the devia- 
tion in entropy rise across the shock wave may be ignored; that is, 

Qvs = Qvs = Qpo or 5 ^vs = 6 Qvs = °- Consequently, an alternate form 

for the attenuation expression 
bp _ Pys ~ Ppo 7 a e 

T 2 a 



where the three terms on the right-hand side of the equation represent, 
respectively, the contributions of region a, of the Q wave generated 
in region p and reflected at the entropy discontinuity, and of the P wave 
generated in region p. Consider now the last term only. The incremental 
change bP along the characteristic is then 
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Sp\ s r tv _1_ 5P 
.^/p vJt f &€ 5t 



(47) 


A change of variable may be made (eq. (2l) of ref. l) to £ since 


oM q - 

U 5t = £ — d| 

oM s - M p - 1 


/sp\ = 1 r l(tv) t^s - Mp 

\ a e/R D if 


/p D J |(tf) ~ M P 


ri c ^,p d * 


'*) - i T Sd - -A - M . e * B c f .(,) At 

^ a €/p D Jo oM s - Mp - 1 ^ 


(48a) 


(&) 

_ 1 


oM s - Mp /0p,pj 

\ a £/(s 

D 

J 0 

oM s - Mp - iy Tp J 


r P c f,P 


(48b) 


Substituting equation (45) into equation (48) yields: 



8p\ 7 a e 1 

V P jp,P " 2 a D 



K(S) r p c f,p 


d| 


(49) 


In equation (49), the left-hand side represents the pressure per- 
turbation at the shock wave due to the wave generation only along the 
forward running (slope of u + a) characteristic. For a complete 
linearized treatment K and Tp may be taken outside the integral and 

equation (35b) may be employed to obtain: 





= Z. — — KTp F 

2 a D PJ 0 


£d 


: f,P 


d£ = Ki 


P ys ~ Pp o) 
Ppo 



P 


(50) 
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Thus, K is the ratio of the contribution of the waves generated 
along the forward running characteristics to the total waves generated 
in region 3 . The remaining portion of the attenuation contribution of 3 
results from § waves reflected at the entropy discontinuity and will be 

designated by (— J A which is equal to 
\ ^ * J3 , Q+S 


, (1 . ( 51 ) 

3,QfS ' p P° /3 



The value of K is plotted against shock pressure ratio in fig- 
ure 7 . The fact that K does not depart significantly from 1.0 means 
physically that the principal contribution to attenuation in region 3 
arises from the P waves. Since figures 5 and 6 have shown region 3 
to have a much larger effect on attenuation in general than region a, 
it is obvious that the theoretical dominating factors for attenuation 
are the P waves of region 3 . This conclusion has been discussed pre- 
viously in references 1 and 2 . 


Since the P waves of region 3 dominate the linearized attenuation 
solution, it is next assumed that a correction for the linearized solution 
may be found by operating only on the P waves in region 3 . Thus P waves 
generated in region a and transmitted at the entropy discontinuity as well 
as the Q waves generated in region 3 and reflected at the entropy dis- 
continuity will retain their original linear or small perturbation values 
even though the attenuation is no longer small. It is further assumed 
that region 3 may be subdivided into a number of smaller regions in each 
of which the linear attenuation relations for P are valid. This treat- 
ment is illustrated in figure 8 . The arbitrary interval A l determines 
the x-wise extent along the shock wave of each of the regions designated 



Each of these regions is bounded by the shock wave 


and two fluid particle paths where each fluid particle velocity is equal, 
respectively, to that generated by the shock wave at the beginning of 
each new interval. The inviscid flow inside each of these regions is 
considered to be constant; and, consequently, there is a small discon- 
tinuity in the inviscid flow across the particle-path boundaries assumed 
in the model. These discontinuities can not, of course, exist in the 
actual physical flow which requires a continuous variation throughout 
all the regions as well as reflections from the shock wave. The errors 
introduced by the assumption of constant quantities in each region are 
not considered to be large and should be of approximately the same order 
as those found in the familiar steady-flow graphical characteristic 
solutions of finite mesh size. 
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In order to simplify the computational procedure, it is assumed 
that for a given region the slope Sx/St of the P characteristic and 
of the shock wave are constant both inside and outside of that region. 

Thus, in figure 8, when the £ contribution of region (T) to attenua- 
tion between A Z and 2 AZ is computed, the assumption is made that 

the characteristic 6,7 and the shock path 0,1 may be extended to. inter- 
sect at point 2. (Numbers refer to points in fig. 8.) The corresponding 


correct regional characteristic line and shock paths are 6,7,8 and 0,1,8 
which are shown in this illustration as intersecting also at the same 
value of x as 2. This intersection at the same value of x is only 
an idealization and is not the true physical picture in general. How- 
ever, since the attenuation effect (generation of £ waves) falls off 
rapidly with distance behind the shock (similar to the fall off in local 
skin friction with distance back of a sharp leading edge in steady flow), 
the contribution to attenuation in the interval from AZ to 2 AZ due 


to generation along 6,7 is much less than that due to generation along 


7,8; thus, small errors in the location of 6,7 will result in very small 
errors in the attenuation at 2 AZ. This assumption for establishing 
the intersection points of the characteristics and the shock wave down- 


stream from a 



without knowledge of the downstream shock-wave 


attenuation permits the easy computation of the influence of region 



for all downstream shock locations. 


When the regional approach described above is applied, the attenua- 
tion for the first interval AZ is identical to the complete linear 
approach. Thereafter, however, the various second-order effects are 

felt. The effect of region (7) on attenuation of the shock during the 

interval from AZ to 2 AZ differs from its effect in the basic linear 
theory because of the convergence of the particle paths since Ujj < Uj. 

This may be shown as follows. From equation (48), 
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The numbers refer to the points on figure 9(a). In this figure the 
lines 5,1, 10,11, and 6,7,2 are drawn with slope (U + a)jj the lines 7,10 
and 0,9,1,11,2 with slope U sI j the lines 0,5,10,6 and 9,7 with slope Ujj 
and the line 1,7 with slope Ujj. Therefore, = 9>7 = 0,10 =. I^O* 


Thus, 




e /i 


1 

D 


r ho 0 M S - Mp 

' loM s - Mp - II 

5 6 \ / 1 


: f0p,i <*1 



(53) 


Thus, the P contribution of region (7) to attenuation between 1 
and 2 is (from eqs. (45) and (53)) 


5p\ 

k p 0/I,ll to l 2 


7 

2 D 






'£) - (54) 

P P/I,0 to 12 \ Pf3 /I,0 to Ijx 


The substitution of equation (40), with g n replaced by g n , into equa- 
tion (54) yields 
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where ( — 


»/l,A Z 


is the linear attenuation of a shock due to P effects 


over the initial interval AZ in region (l)- 


The relationship between Ipp and is derived in appendix B 

and is 


2 11 
2 1 



Consequently, 


8p\ 

P °°/I, Z^ to Z 2 




&p\ 

Po °/I,Z 2 to Z3 



n+1 
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(56) 
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which for equal intervals AZ becomes 


*\ 


5 p 


_ / 6 P 


■»/i } l ± to Z 2 \ °°/I,AZ 


' 5 p 
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I,Z 2 to Zj \ °°/I,AZ _ 
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n+1 
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n+1 

Zij n+ 3 


- H 


n+1 

n+3 


(57b) 


In this form the nonlinear attenuations can easily be computed. 

In this manner the influence of region ( 1 ) on attenuation at any 

desired value of Z may be computed once the attenuation at Z^_ has 
been found. The influence of region (il) is found in a similar manner 

by shifting the effective origin of the coordinate system to Z^ and 

finding (— ) for the attenuated shock strength at Z-|_; that is, 
\ Po °/II 


(Sp 


( 5p\ 

^P- 


XX y Zj^ "to z 2 
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(P„ 


'II,AZ 


II, Z 2 to Z3 \ P “/II,AZ|_' 
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(58a) 
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or for equal intervals of A 2 , 


Sp 


6p 

Poo } 


■\ 


III, l 2 to Z 5 \ “y II,AZ 


n +1 

2 n+3 


" Pi! 


n +1 

n+3 


f =f 
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n +1 

jn+3 


-( 


2Z II 
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( 58 c) 


The total nonlinear attenuation at a distance l from the diaphragm 
station which has been subdivided into several intervals A2 is then 
expressed as the sum of the linear contributions of region a and the 
reflections from region P added to the nonlinear contributions of region p. 
The following expression is obtained for the nonlinear attenuation: 



where i is the index. 


(59a) 


(59b) 
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Second-order attenuation in the presence of transition may also he 
treated by this regional system. The contributions of a and of Q 
and S in P are still treated as entirely linear but including transi- 
tion. Equation (59a) applies in its entirety to the a contribution; 
and, when the first term on the right-hand side is multiplied by 1 - K m 
and the second term by 1 - K n , the resulting equation gives the Q 
and S contribution of region p. 

The transitional P contribution is treated as follows: Let Z* 

be the position of the shock at which the transition in the flow behind 
the shock first affects the shock. The wave diagram for transition is 
shown in figure 9(b) where for illustrative purposes it is assumed that 
Z* = 2 AZ. Now, in order to retain the facility of computation afforded 
by the regional system with equal AZ and a constant Z*, the value of 
R* will, as a result, vary slightly from region to region. Since 

the magnitude of this variation may be found by examina- 
tion of figure 10 which shows the parameter 



U p oo, std 

V Z*/ N Poo 


plotted against shock pressure ratio. From this figure it is evident 
that, if the shock pressure ratio should attenuate, for example, from 20 

to 15, or from 5 to 4^, there would be about a 10-percent decrease in 

2 

R* for a given Z* and p^ std ^P 00 * 1116 errors introduced by such a 

variation in R* are not deemed to be important enough to force the 
abandonment of the equal AZ computing scheme. For the remainder of 
this paper, R* will be taken as the value of the transitional Reynolds 

number in region The Sp contributions for the transitional case 

are expressed by the following equations which are modifications of equa- 
tions (57), (58), and (59) (the subscripts m and n refer to condi- 
tions before and after transition, respectively): 


S£ 

l Poo /N,m,AZ 



m+1 

m+3 


(60a) 


36 


NACA TN 4347 



N,n,AZ 



2 n+1 



( 60 b) 
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(6lb) 


for 



N 


1 + 


if 

AZ 


Numerical Evaluation of Nonlinear Theory 

Several computations to determine the nonlinear correction factors 
for attenuation vere performed for values of the interval AZ/D of 3*5 
and 14. These particular values vere chosen because they represent 
increments of AZ of O .5 foot and 2 feet for the shock tube employed in 
the experiments to be described in a later section. Typical curves 
resulting from such computations are shown in figure 11. The ratio of 
the nonlinear attenuation to the linear attenuation is plotted against Z/D 
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for an initial shock pressure ratio of 4 and a value of 
For values of Z/D > 50 the curves for ^ = 3«5 and. l4 


of 0.01 X 10 6 . 

v co 

agree within 


about 1 percent for the laminar case and within about 3 percent for the 
turbulent case. The low value of a^D/v^ accentuates any variations 
between the two computations; thus the case illustrated gives a discrep- 
ancy near a maximum rather than near a minimum. Examination of several 
such pairs of curves resulted in the conclusion that the slight ^increase 


in accuracy obtained by using — = 3*5 did not justify the fourfold 

increase in labor. Consequently, the computations with — = l4 are 

D 

used to predict the nonlinear attenuation for z/D 5 0 . 


It is obvious that the finite size of AZ will introduce errors 
in the ratios NIm,n/lm,n which are largest near Z -* 0 since the non- 
linear and linear attenuations are identical for the first interval. 

(See fig. 11.) However, the errors introduced in the attenuations Nl^n 
themselves are small since I^n -» 0 as Z -»0. To represent the physi- 
cal flow in this region accurately would require that AZ approach 0. 
Interpolation formulas giving acceptable accuracy near Z -*0 are assumed 
to have the form 


Nonlinear attenuation 
Linear attenuation 


(Constant )^- 


0.5 


Nonlinear attenuation 
Linear attenuation 


= (Constant)^ 







(62) 


for ^0 < ^ ^ 56 j for laminar and turbulent flows, respectively, since 

the linear attenuation is proportional to (z/d)^'^ and ( Z/D) for 

the laminar and turbulent flows . 

Al 1 1 

The constants are chosen to match the computed curves for — = 14 


at i = 56 . From figure 11 it is evident that the errors resulting from 
the application of this interpolation formula are less than the afore- 


mentioned errors at 



and. are thus acceptable. 
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In order to obtain curves of the transitional nonlinear to linear 
attenuation ratio NL^t/I^T f° r a constant value of R* and several 

shock pressure ratios, a cross -plotting technique was used. At a given 
shock pressure ratio the values of NL7 ,T were computed for several 

values of the ratio Z*/D to AZ/D. Since each Z* represented an R*, 
the resulting ratios for each Z/D could be plotted against R* for a 
given shock pressure ratio. The values for a particular R* could then 
be read from these plots to produce a master plot with a common value 
of R*. 

Plots of the ratio of nonlinear attenuation to linear attenuation 
are shown in figures 12 to 15- In these curves region a is always con- 
sidered as having turbulent flow whereas four cases are considered for 
region p: namely, (a) laminar flow, (b) turbulent flow, (c) transition 
with R* = 1.25 x 10°, and (d) transition with R* = 2.5 X 1C>6. The 
cross -plotting parameters are shock-tube Reynolds number a<oD/ and 
length of shock-wave travel expressed in hydraulic diameters Z/D. At 
the lower pressure ratios Ppo/Pco> curves for more values of aoJl/voo are 

shown them at the higher pressure ratios. This limitation resulted from 
the considerations of the restriction of the validity of the theory to 
an ideal gas, the region of experimental data of this report, the most 
likely generel region of experiments for other facilities, and priority 
for computing effort. Since figures 12 to 15 are the result of cross 
plotting, the accuracy is assumed to be about 2 percent. 

An analytic closed-form investigation has been made of the fact 
that the limit of approaches 0 when Pp 0 ^/P ro approaches 1 

whereas the limit of Hl/j jljj approaches 1.0 when Pp 0 ^P ()0 

approaches 1.0. This second-order analytic solution for weak shocks 
indicates that a value of n = 3 in P is a critical value; all solu- 
tions with n > 3 approach a limit of 1 and those with n < 3 approach 
a limit of 0. Of course, in all cases the absolute value of both the 
linear and nonlinear attenuation must approach zero as Vfi 0 jp m 

approaches 1. Since a has only a secondary effect on attenuation and 
since the expansion fan has been replaced by a "negative shock," the 
refinement of transition in region a was not deemed necessary. 


RESULTS AND DISCUSSION 


Experimental Apparatus and Procedure 


Extensive shock-attenuation data were 


shock tube 2 inches high by 1^ inches wide 


obtained in a high-pressure 
in the Langley gas dynamics 
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laboratory. This shock, tube is the same as that described in reference 1 
with added velocity-measuring equipment. Light screens placed at eight 
stations in the low-pressure section made possible measurements of the 
complete distance -time history of the motion of the shock wave for a wide 
range of shock pressure ratios and flow Reynolds numbers. Figure 1 6 shows 
schematically the arrangement of the shock tube and associated equipment, 
and the low-pressure section with associated optical systems is illustrated 
in figure 17* 

Air at room temperature was used in both high- and low-pressure sec- 
tions. The normal arrangement for high-pressure section air supply and 
low-pressure section vacuum systems is shown in figure 1 6. In a limited 
number of low-density runs, evacuating the high-pressure section was nec- 
essary; for these cases, an auxiliary vacuum system, identical with the 
normal low-pressure system, was substituted for the pressure system shown. 
Conversely, certain high-density runs required pressurization of the low- 
pressure section; in this case a simple pressure system replaced the nor- 
mal vacuum system. The pressures in both high- and low-pressure sections 
were adjusted for each run. All data were obtained from tests where the 
diaphragm was punctured by a hand-operated plunger. In this way diaphragm 
pressure ratio and consequently, theoretical shock-pressure ratio Pp 0 /Poo 

were closely controlled. Bourdon-tube gages were employed for pressure 
measurements, and the vacuum systems utilized a modified barometer for 
pressures in the range from 4 inches mercury absolute to atmospheric 
pressure, and a 0 to 100 millimeters mercury absolute pressure gage was 
used for the low pressures. 

Static pressures were structurally limited to 1,000 pounds per square 
inch gage in the high-pressure sections and 250 pounds per square inch 
gage in the low-pressure sections. A vacuum limitation of about 0.01 atmos- 
phere absolute pressure resulted for the low-pressure section because the 
light screen systems became unresponsive for the low pressures. In the 
high-pressure section, low pressures were limited by failure of diaphragms 
to burst properly. 

Diaphragms made of thin metal foil were used for the low-pressure runs. 
The most useful materials were soft aluminum foil, 0.001 inch thick, and 
soft brass foil with a nominal thickness of either 0.00125 or 0.0015 inch. 
With these materials, it was possible to obtain good bursts for pressure 
differences across the diaphragm ranging from 10 pounds per square inch to 
100 pounds per square inch. Good bursts for the range of pressure differ- 
ence from 100 to 1,000 pounds per square inch were obtained by using spring- 
tempered brass shim stock with thicknesses ranging from 0.008 to 0.021 inch 
and scribed to various depths in an x-shaped pattern along the diagonals of 
a rectangle representing the shock-tube cross section. When punctured 
under pressure, the diaphragm split along the scribe marks, and the four 
triangular pieces of material folded back against the wall and presented 
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minimum resistance to the flow. All other conditions being equal, r\ons 
where this type of burst took place resulted in minimum shock attenua- 
tion compared with attenuations resulting from bursts where metal or 
acetate-type diaphragms were shattered. Material thicknesses and scribe 
depths were determined so that puncturing pressure was just below the 
pressure at which the diaphragm would have failed. The unscribed foil 
diaphragms used for small pressure differences split along the diagonals 
in this same way when punctured at the center, and optimum results were 
also generally obtained with these diaphragms. 

Shock velocity measurements were made with miniature schlieren systems 
located at eight positions in the low-pressure section. Figure 1 6 shows 
schematically the position of these systems. The optical and electronic 
systems were essentially the same as those used in reference 1, in which 
the signal generated in a photomultiplier tube by deflection of a beam of 
light upon the tube was amplified and used to trigger a thyratron. The 
thyratron output pulse, in turn, started or stopped a counter chronograph. 
Figure l8 illustrates one complete optical system, including the chassis 
containing the photomult iplier-amplif ier-thyratron circuit, which is shown 
in figure 19. Wooden shields were employed to keep stray room light from 
falling upon the photomultiplier tube. 

For the multiple systems employed herein, each thyratron output pulse 
was channeled to two chronographs. These chronographs indicated shock 
traversal time between any two adjacent stations, and any one of three 
types of chronographs, 8 megacycles, 1.6 megacycles, and 1 megacycle, 
was employed for each position. The overall time interval between the 
first and last stations was measured with a 100-kilocycle chronograph 
for a check upon the sum of the individual measurements. 


The experimental shock-wave pressure ratio was computed from the 
measured time interval At for the shock to pass between two measuring 
stations a distance AZ apart from the relation 


This value was assumed to represent the shock strength at a position 
midway between the two stations, and the maximum error was estimated to 
be less than 1 percent. The theoretical shock pressure ratio Pp 0 /£oo 

was computed from the diaphragm pressure ratio just prior to burst. For 
an ideal gas the maximum error in pp Q was estimated at O.lp^ ^for a 


Reduction of Experimental Data 
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shock pressure ratio of 10 at = O.CKTpp^ gtd J; 8X1(1 the maximum devia- 


tion in the ratio of computed 
to be approximately 1 percent. 


( P p/ P “)o 


to true 


N p »)o 


was estimated 


Comparison of Theory and Experiment 

The general method of comparison between theory and experiment will 
be to compare the measured and predicted attenuation for particular values 
of Ppo/Poo and- a <xP/ v <x> on individual curves. However, it is of interest 

to first consider a few typical curves where the data for a constant value 
of Pp^p^ but with varying values of a 0 fi/ v C o 8X6 shown on a single plot. 

Such plots are shown in figure 20 for values, of Pp 0 /P M of 4.0 and 10.0. 

The experimental data are averages of several runs on a given day and the 
number of runs for each data point is indicated on the figure. On some 
runs in which the density change across the shock wave was small, all 
the velocity-measuring stations did not register because of variations 
in their sensitivity; and, as a result, there are gaps in the experimental 

data, (fot example, see fig. 20(b) at = 0 . 005 ^ 

One fact immediately evident is the nonrepeatability of some of the 
data, even when compared on a daily average basis. An allowance for an 
error in p^p^ of ±1 percent in the experimental and data-reduction 

technique will bracket most of the observed discrepancy in the averages; 


but certain runs at very low values of 
range. 


sWD f\ 


still fall outside this 


If an attempt were made to extropolate a curve from the data points 
to l = 0, an inflection would often be required in the curve between 
1=0 and 1=6 to make it pass through the theoretical value of 
P(3 o/ Poo at 1 = 0. In order to illustrate this point, connecting lines 

have been drawn in figure 20 for some of the values of a^D |v M . Similar 

behavior is found in the experimental data reported in figures 11 to 14 
of reference 8. Now all the attenuation theories based on wall effects 

^Pvs/Poo 


which are known to the authors predict 


>0 in regions of laminar 


dV 


or turbulent flow. At the transition point the theory of the present 






< 0 . Con- 


paper usually predicts — — x — ) > . 

'laminar \ /transition 

sequently, if this inflection is to arise from wall effects, it must be 
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caused by transition. The interferometric studies of reference 16 have 
found values of R* of about 1.4 x 10^ and 2.0 X 10^ for Pp^P*, = 2 

and- 24, respectively. Since values of R* of the order of 0.5 X 10^ 
are required to cause an inflection at Z = 6 for ^ and 10, 

it does not appear likely that transition is the cause of this inflec- 
tion. This inflection will be discussed more fully later in this paper. 

The unexpected inversion of Reynolds number effects for the lower 
values of Z should be noted. For example, at Z = and Pp^P^ = 10 

the attenuation for = 0.005 x 10^ is about one-fourth that for 

Voo 

= 0.1 x 10^. For the same value of Z and ^2 = 4.0, the attenua- 
v » Poo 

tion for = 0.1 X 10^ is about one half that for = 1 x 10^. 

v„. 


The spread with Reynolds number of the experiments! data at the 
larger values of Z is also much smaller than would be expected on the 
basis of the linear theories (refs. 1 and 2) which predict attenuations 


at a given Z/D proportional to 




1 

5 

for laminar 


and turbulent flow, respectively. This behavior has been noted by other 
investigators (ref. 10). 


Theoretical variation of p, /p with Z for values of 

r vs/ ^00 

= 0.01 X lO^, 0.1 X 106, and 1 X 10^ for p R -/p = 4.0 and 10.0 are 
Voo P °/ 00 

presented in figures 21 and 22. These values are based on the theory 
of the present paper. Laminar and turbulent linear and nonlinear 
pressure-distance predictions are shown in figure 2l(a) for Pp^P^ = 4.0, 

and in figure 2l(b) for = 10.0. The effect at larger values of Z 

of the nonlinear correction is twofold; not only is there a marked reduc- 
tion in the predicted attenuation when the attenuation is large but 
there is also a marked decrease in the predicted Reynolds number effect 
on attenuation. The following table based on figures 21 and 22 for 
Z = 22 feet illustrates these facts: 
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Curve 

Figure 

Values of for ^ = 4 

a_D 

for values of - — of - 

Values of for ^-10 

\ P «/ \ P ® / Z=22 P ® 

a»D 

for values of - 2 - of - 

R* 

0.01 X 10 6 

0.1 x 10 6 

lxiO 6 

0.01 x 106 

0.1 X 10 6 

1 X 10 6 

*7,1 

21 

1.45 

0.48 

0.17 

2.62 

0.66 

0.10 

00 

*7,7 

21 

2.30 

1.45 

• 91 

6.10 

3.85 

2.43 

0 

NL7,1 

21 

1.09 

.44 

.16 

1.87 

• 58 

.10 

00 

nl 7>7 

21 

1.39 

1.04 

.74 

3.22 

2.50 

1.83 

0 


22 

1.09 

•99 

• 73 

2.65 

2.46 

1.8 

1.25 X 10 6 


22 

1.09 

.92 

.72 

1.87 

2.38 

1.8 

2.5 X 108 


The nonlinear attenuation predictions for values of R = 1.25 X 10° 
and 2.5 X 10^ are shown in figures 22(a) and 22(h). Since, for these 

values of R*, Z* = 0.122 and 0.244 foot for = 106 at ^§2 = 10, 

V 00 Poo 

the nonlinear curve for turbulent flow is nearly identical to that for 
transition and the turbulent curve is used in the figure. These curves 
and the table indicate that transition at a constant value of R appears 
to decrease further the spread of attenuation with The discon- 

tinuous change in d Py . s / P ” at Z* is obvious. Another point of 

dZ 

interest is the predicted variation of NIz^^ at Z = 22, Pp 0 /Poo = 10* 

and R* = 2.5 x lO^j the attenuation is less for a^D/v^ of 0.01 x 10^ 
and 1 x 106 than it is for a^/v^ of 0.1 x 10 ^. The completely laminar 

flow of = 0.01 X 106 and the nearly completely turbulent flow of 

Vqo 

= 1 x 10^ result in almost identical attenuations. The flow of 

Voo 

fsB = 0.1 x 10^ is mainly turbulent ( Z* = 2.44) and, as expected, the 

Voo 

attenuation is greater than for = 1 x 10 ^. Thus, in this case, 

V 00 

transition has resulted in an inversion of attenuation with Reynolds 

v 6 

number under a certain set of conditions. However, for R = 1.25 X 10°, 

= 10 , and Z = 22 feet, the pattern of increasing attenuation with 
Poo 

decreasing a^D/vo, is once again evident. 
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In figures 23 to 36 the variation of the average experimental shock 
pressure ratio with distance is compared with the various attenuation pre- 
dictions of the theory presented earlier in this paper. The nominal 
values of P(3 0 /Poo range from 1.5 to 10 and a^Djv^ covers the range 

of 0.005 x lO^ to 15 X lO^. The hydraulic diameter of the shock tube 
is l/7 foot and the maximum value of Z at which Vys/^oo cou ld be 

Q 7 

determined was feet; thus a maximum experimental value of ^ ** 125 

resulted. 


The marked improvement obtained from use of the nonlinear theory 
when the predicted attenuation is large becomes evident upon inspection 
of figures 25 to 56. The unusual behavior near Z -> 0 which was men- 
tioned earlier can now be examined more closely. For the higher values 

ra T) S' 

of shock pressure ratio (8, 9 , and 10) and low values of < 0.01 x 10 , 

^00 

the measured attenuations at the first station are much smaller than any 
of the theoretical predictions of this report or of reference 2. In fact, 


for ^ = 9.0 
Poo 


and 


— = 0.005 X 10 6 , 
^00 


one set of averages gives a nega- 


tive attenuation D f such size that even the estimated 

\P« Poo / 


Pvs 


1-percent-error margin is not sufficient to make — — < 

Poo Poo 


po 


at the first 


( \ CO -^00 

^ = 56. 6j. In reference 9 if was also found that for high values 
of Pp 0 ^Poo the maximum shock velocities were greater than the velocity 


theoretically computed for an inviscid fluid both with and without con- 
sideration of variable specific heat and gaseous imperfections. Hydrogen 
and helium were used as the driver gas and air as the low-pressure fluid. 
The maximum shock velocity occurred at about 40 to 50 diameters from the 
diaphragm. 


On the other hand, references 4 and 8 did not find any values of 


> Ej3o f or values of 
P ro Poo 



in the same range as the present 


experiments. These works proposed a "formation decrement" defined as 
"the difference between the Rankine-Hugoniot shock strength and the 
maximum shock strength obtained after the formation distance." (See 
ref. 8, page 17*) This decrement was then attributed to the imperfect 
diaphragm burst producing a series of compression waves which eventually 
coalesce to form a shock weaker than that for the case of a theoretically 
perfect burst. 
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For the lowest shock pressure ratios at the highest shock tube 
Reynolds numbers, the opposite trend (measured attenuation larger than 
theory) is apparent. An extreme example of this is the attenuation meas- 
ured at the first station for = 15 X 10^ and = 2. (See fig. 25 .) 

V «. Pco 

Three of the four runs used in obtaining this average value showed this 
behavior, which is attributed to the formation process. 

The processes giving rise to the behavior for Z/D -»0 arfe not 
understood. The experimental data of this report generally show a trend 
of decreasing P vs /P|3 0 with increasing a^D/v^, at the first station. 

If this effect were to be explained on a viscous basis, it would appear 
that an opposite trend would appear; hence, the diaphragm burst is prob- 
ably the governing factor. The bursting phenomena are in turn governed 
by the diaphragm materials (weights and rupture characteristics) as well 
as the pressure load. The diaphragm opening time (time for the diaphragm 
sectors to fold against the wall) was estimated by a method which has 
previously been checked experimentally. No correlation was found between 
the opening time and the behavior near Z — »0. It was noticed, however, 
that, in the experimental runs which exhibited the most marked inflections, 
foil diaphragms were used. On the other hand, some runs exhibiting very 
little or no inflection also used foil diaphragms. Therefore, it does 
not appear that foil diaphragms are solely responsible for the inflected 
data points. 


Regardless of whether the maximum experimental value of Pys/P,,, I s 

greater than or less than the ideal value, the behavior near the diaphragm 
station is not governed by wall effects but by the diaphragm burst and 
the resultant three-dimensional flow first established. This highly rota- 
tional viscous flow does eventually become essentially two-dimensional, 
with the exception of mixing and vorticity at the interface between the 
driver and driver gases. Consequently, any attenuation theory based on 
wall effects cannot predict the initial behavior near the diaphragm sta- 
tion. As the distance from the diaphragm station increases, the ratio 
of the influence of the initial bursting flow to the influence of the 
integrated wall effects decreases; therefore, the physical variation of 
Pvs/Poo should approach the theoretically predicted attenuated value 

asymptotically as z/D increases. 

The experimental data of figures 23 to 36 approach the nonlinear 
theory (considering transition) in a manner very similar to that just 
described. For high values of Vp 0 /V m and low values of a^D/vo, where 

the effect of the formation process results in a significant inflection, 
the deviation from the nonlinear curve persists to the larger values of 
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Z/D. For — i — = 10, 9, and 8 with = 0.005 X 10°, the formation 

Poo V oo 

effect was so large that the experimental data never approached closely 

the nonlinear curve for ^ ^ 125* As a^jv^ increased and the forma- 
tion effect decreased, agreement between experiment and nonlinear theory 
improved both in convergence of the experimental and theoretical values 
at a lower value of Z/D and in maximum deviation at the highest ' values 
of Z/D. 

The nonlinear curve for R* = 2.5 X 10^ appears to agree more 
favorably with the data over most of the range of Ppo^P^,- At the lower 

values of P^q^P^ there appears to be a tendency for the data to depart 

from the R* = 2.5 X 10^ curve and approach the R* = 1.25 x 10^ curve. 
The interferometric measurements of reference l6 indicate such a trend 
of increasing R* with increasing Ppo/Poo’ This trend might also be 

expected from comparison with steady-flow experiments since the wall 
cooling increases as Ppo^P^ increases. 

The comparisons of figures 23 to 36 between the nonlinear transitional 
theory and the experimental data show that this theory is valid for the 
prediction of experimental attenuation except for the lower shock tube 
Reynolds numbers at the higher pressure ratios. These latter conditions 
are those under which it appears that the shock-formation processes dom- 
inate the entire flow. 

As stated previously, errors of unknown magnitude were introduced 
by the averaging of wall effects across the flow. The fact that the 
present theory was able to predict fairly well the measured attenuations 
over the entire range of shock pressure ratios and Reynolds numbers 
(except near the diaphragm) indicates that the errors introduced in the 
averaging process are either not serious or else self- compensating. 

The range of boundary-layer thicknesses in region (3 (which has the 
predominant influence on attenuation) was determined by methods similar 
to those of references 14 and 15 for laminar and turbulent flows, respec- 
tively. In the experiments reported herein, with the shock 20 feet from 
the diaphragm, the maximum boundary-layer thickness varied approximately 
from 0.006 to 0.25 inch for laminar flow and from 0.07 to 1 inch for tur- 
bulent flow. When compared with the 0. 75-inch half -width of the shock 
tube, the turbulent thickness was significant for much of the experimental 
data. On the other hand, laminar boundary-layer thickness was generally 
small. Thus the averaging process as used for this theory appears to 
apply equally well to thick as well as to thin boundary layers. 
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Magnitude of Other Neglected Second-Order Effects 

Since the treatment presented has considered only one source of 
nonlinearity, it is advisable to examine briefly the influence of other 
neglected sources. The three most obvious factors neglected are the 
finite extent of the expansion fan, the nonlinear effect of the wave 
interactions at the entropy discontinuity, and the reflections at the 
shock wave. 

The effect upon attenuation of treating the expansion wave as one 
of finite width with varying free-stream properties was calculated by 
using the results of reference l4 for pressure ratios Vp 0 /l>oo of 1 * 6 > 

2 . 9 , and 4.5 (tte of 8111(1 respectively, in the notation of 

ref. l4). Computations were not made for higher shock pressure ratios 
because the solutions in reference l4 did not extend above | te = 0 . 6 . 

The laminar-flow wall-shear and heat-transfer distributions through the 
cold-gas regions given by this reference were used to compute the skin- 
friction coefficient which was, in turn, utilized to compute the 
attenuation. 

When the method of the present paper was used, it was found that, 
for the three cases computed, the net change in attenuation through 
approximation of the finite expansion by the negative shock was practi- 
cally zero for Pp 0 ( /P OT equal to 1 . 6 , - 0.3 percent for P^P,*, equal 

to 2 . 9 , and -0.6 percent for P^o/p,*, equal to 4.5. Although the con- 
tribution of the cold gas was itself influenced by the finite expansion 
(up to a 50 -percent decrease at the highest pressure ratio), the cold 
gas contributes so little to laminar-flow attenuation at these low pres- 
sure ratios that the error in assuming a negative shock is negligible. 

For higher shock pressure ratios the effect of the finite expansion 
cannot be computed. However, examination of figure 5 indicates that, at 
a shock pressure ratio approximately equal to 6 , the contribution of the 
cold gas (with the negative shock) vanishes. At higher pressure ratios 
the cold gas tends to decrease attenuation but this effect remains small 
compared with the contribution of the hot gas. At a shock pressure ratio 
of 20 , for example, the effect of the cold-gas region has reached only 
4 percent of the total. Thus, the finite expansion fan can influence 
only a gmai 1 part of the total attenuation, and the assumption of a nega- 
tive shock should give reasonably accurate results. 

The laminar boundary-layer finite-expansion-fan solutions of ref- 
erence l 4 were also used to estimate the effect of the negative shock 
assumption upon the attenuation predicted by reference 2 for the same 
three pressure ratios. The errors in attenuation which arise through 
the use of a negative shock are 1.4 percent at Ppo/^oo e 9 ua d to 1 . 6 , 
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5.1 percent at V^q/v^ of 2.9, and 4.2 percent at p^p^ of 4.5. In 

this case, the contribution of the cold gas tends to decrease the attenua- 
tion at all pressure ratios of interest, the contribution for the negative 
shock assumption reaching a maximum at Pp 0 /Poo e H ual to 

No means of computing turbulent boundary- layer flow inside a finite 
expansion are available. However, according to the linearized negative 
shock approach of this paper, the cold-gas contribution is always less 
than l6 percent for pressure ratios up to 20 and always less than 10 per- 
cent for pressure ratios up to 10. (See fig. 5>) If the effect of the 
finite expansion in turbulent flow is in the same direction as in laminar 
flow, the cold-gas contribution, when a finite expansion fan is considered, 
would be less than that for a negative shock. Thus it is assumed that the 
effects of the difference between the finite -expans ion and negative- shock 
solutions for turbulent flow may be neglected., at least up to shock pres- 
sure ratios of 10. 

The influence of nonlinearities in the reflection and transmission 
of waves at the entropy discontinuity was calculated for shock pressure 
ratios Pg 0 |Poo equal to 1.25, 2.0, 4.0, 6.0, 8.0, and 10.0. The non- 

linearities involved required that the shock-tube Reynolds number and 
station be specified in order to calculate the attenuation for this case. 
The case of equal to 0.1 X K)6 was chosen, and the attenuations 

at Z/D equal to 70 and 140 were calculated for completely laminar (n = l) 
flow and completely turbulent (n = 7) flow. 

For the laminar case, a difference of. about one-half percent was 
found between linear and nonlinear calculations for the wave interactions 
at the entropy discontinuity for of 10 and at Z/D equal to l40, 

the nonlinear calculation predicting greater attenuation. Differences at 
Z/D of 70 and at lower pressure ratios were smaller and generally in the 
same direction, except at a value of Ppo/Poo of 1*25 where the nonlinear 

entropy result gave slightly less attenuation than the linear result. 

The deviations, which were small in all cases, were also somewhat erratic 
in their behavior. 

Differences in the turbulent case were somewhat larger and ranged 
up to 10 percent for a value of Z/D of l40 and a value of Ppo^Poo 0 I 10* 

This condition is attributed to the larger relative contribution of the 
cold-gas region a to the total attenuation for turbulent flow. (See 
table I and fig. 5.) Again, the nonlinear entropy calculation predicted 
greater attenuation than the linearized except at a value of Pp G /Pco 

of 1.25 where the trend was reversed. However, at the low pressure ratios 
the deviation was less than 1 percent and this trend is not considered to 
be significant. 
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The effect of reflection of waves at the shock was estimated by the 
following treatment: The unattenuated shock was assumed to be overtaken 

by a single isentropic wave of strength equal to the total attenuation 
for a given condition and station. The resulting one-dimensional inter- 
action was computed and the strength of the shock after interaction com- 
pared with the strength predicted by the linearized theory. These results 
were computed for shock pressure ratios of 10 and 4. 

For Ppo^Poo equal to 10, the three cases were compared and the 
results Eire given in the following table: 



z 

D 

Attenuation 

Percent difference 
in attenuation 

Voo 

No reflection 

With reflection 

o.i x io 6 

140 

1^ 1 = - 0.60 

L i,i = - 0.55 

8 

.1 

140 

L l, 7 = " 3 * 57 

L 7^7 = - 3 • 16 

11 

lf\ 

O 

O 

154 

L7,7 = _7 -° 

lyj rj = -6 . l4 

12 


The conditions of -4 — = 0.005 X 10^ and — = 154 represent the condi - 

D 

CO 

tions under which maximum attenuation would he expected for the experi- 
mental range of this paper. 


In all cases, the shock is slightly strengthened by the interaction, 
that is, attenuation is decreased. This is a trend which is opposite to 
that generally computed for the exact entropy discontinuity. For Pp 0 /Poo 

equal to 4, only one case was computed, that corresponding to the largest 
predicted attenuation, a^J) equal to 0.02 X 10°, l/D equal to 154. 

In this case, the predicted attenuation was ^ vs ~ equal to -2. The 

Poo 

consideration of reflection gave vs — ^ equal to -1.8. The differ- 

Pco 

ence in attenuation was 10 percent, again in the direction of decreasing 
attenuation. 

In general, then, the effects of the entropy discontinuity and the 
reflected wave at the shock which were neglected in the linearized theory 
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tend to oppose one another and are of about the same order of magnitude 
for the worst cases ( turbulent flow, strong shocks). 

For the pressure range of this report it has been shown in the pre- 
ceding discussion that the two nonlinear effects having the largest magni- 
tude are of opposite .sign. The effect of neglecting the finite expansion- 
fan width is small for laminar flow and is also presumed to be small for 
turbulent flow. Consequently, the neglect of these three effects, which 
have a net influence much smaller than the nonlinear P effect considered, 
appears to be justified for attenuation. Although the aforementioned 
effects are negligible for attenuation, they can appreciably influence 
the pressure, density, and velocity distribution at points removed from 
the shock. This behavior arises from the fact that the relative influence 
of the £ waves of region |3 decreases as the distance behind the shock 
increases; the influence reduces to zero at and behind the entropy dis- 
continuity. (See pages 35 £uid 3^ and fig. 13 of ref. 1 for further 
discussion. ) 


CONCLUDING REMARKS 


The theory of NACA Technical Note 3575 , in which shock-wave attenua- 
tion is calculated by use Of a linearized form of the hyperbolic system 
of equations of motion and energy and through the assumption of equivalent 
incompressible steady-flow friction and heat-transfer coefficients, has 
been modified in the following manner: 

(1) Incompressible unsteady skin-friction coefficients have been 
determined by an integral method. The resulting unsteady incompressible 
skin-friction coefficient is corrected for compressibility by a refer- 
ence temperature method. 

(2) A nonlinear regional approach has been employed to permit the 
extensions of the theory to large attenuations. This approach modifies 
only the forward running waves generated in the hot-gas region; these 
waves are shown to dominate the attenuation process. 

(3) Transition effects are considered. The method assumes instan- 
taneous transition from laminar to turbulent flow. The Reynolds number 
of transition then becomes a parameter of the attenuation problem. 

The modified theory has been evaluated for a range of shock pres- 
sure ratios of general experimental interest. Curves are presented to 
permit easy prediction of attenuation for shock pressure ratios to 20 
and a range of shock -tube Reynolds numbers. Results for the linearized 
theory with all laminar and all turbulent flow are compared with results 
of NACA Technical Note 3278 , and the predicted attenuations are found to 
be in fair agreement. 
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Experimental local shock pressure ratios have been determined for 
ideal shock pressure ratios from 1.5 to 10 for a range of shock-tube 
Reynolds numbers and position from the diaphragm station. 

Comparison of the modified nonlinear theory using a transition 
Reynolds number based on flow length of 2.5 X 10° with the experimental 
results shows good agreement except for the following situations: 

(1) The highest shock pressure ratios at the lowest shock- tube 
Reynolds numbers where the effects of the nonperfect diaphragm burst 
are believed to dominate the flow. 

(2) The lowest shock pressure ratios, for which a lower transition 
Reynolds number (l.25 x 10o) appears to give better agreement. 

The effects of considering a finite expansion fan, the exact inter- 
action at the entropy discontinuity, and reflection at the shock wave 
(all of which are neglected in the present treatment) are calculated for 
certain cases. The net effect of these three contributions is shown to 
be small compared with the nonlinearities accounted for in the stepwise 
regional calculation. 


Langley Aeronautical Laboratory, 

National Advisory Committee for Aeronautics, 
Langley Field, Va. , June 25, 195®* 
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APPENDIX A 


SUMMARY OF LINEARIZED ATTENUATION EQUATIONS 

The changes in the characteristic parameters P, Q, and S/R as 
a function of local skin friction are derived in reference 1 and those 
appropriate to shock -wave attenuation are summarized in this appendix. 
Subscripts refer to positions on the wave diagram of figure 1. 


p q _ = . p c ~ Pq ° = 0p a -^4 _ + T] — 1 r ld c fa (£)d| (Al) 
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APPENDIX B 

DERIVATION OF EQUATION FOR % 


The derivation for Zq will he obtained and then generalized to 
arbitrary N. The wave diagram of figure 9(a) shows the location of the 
numbered points. The symbol xg denotes the distance of the point 6 

from the origin, and so forth. The following equations are obtained by 
s imp le geometry from figure 9( a ): 


x 2 = U sI t2 

x 2 = U s it]_ + Uxx(ty - tx) + (U + a)x(t2 - ty) 


(Bl) 

(B2) 


Combining equations (Bl) and (B2) yields 


^7 = 

L 1 u sl - U II 

^2 _ 

u sl ' U II 

ti 

Uji - (U + a)j 

t l 

Uji - (U + a)j 


Also from figure 9( a )> 


(B3) 


X 11 = U sl^ll (B4) 

Xu = U sI t 1 + U-^ty - tx) - U sI (ty - tio) + (U + a) I (t 11 - tqo) (®5) 

x n = u ltio + (U + a ) I ( t ll - t io) ( b6) 


Combining equation (Bl) with equation (B6) and combining equation (B5) 
with equation (B6) yields 


(B7) 
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and 

tio = u si - "a h _ A (b8) 

tl u sl " U I \Jl / 


Now 


h± _ ^11 = ^11 ^10 

Z 1 t l t 10 t l 


(B9) 


Substitution of equations (B3), (B7)* and. (B8) into equation (B 9 ) pro- 
duces 


2 11 = U sl ~ U II a I A 2 _ A 

2 1 U sl ~ U I (U + ajj - Un^t! J 


or 


Z 11 

2 1 



(BIO) 


(Bll) 


Since the origin in figure 9 could be shifted to any arbitrary 
location and the same geometrical relations above could be derived in 
the new location the generalized formula for Zjf can be written from 
inspection of equations (BIO) and (Bll) as 

„ u sN " U N+1 a N 

Zjj 

u sN " % (U + a) N - %+i 


(B12) 
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Figure 1.- Distance-time plot of shock-tube flow. 
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Incompressible laminar skin friction, Cf 



Figure 2.- Plot of incompressible local skin-friction coefficient for wave-induced, flows having 
ratios of wave velocity to flow velocity corresponding to those in a shock tube of given 
shock pressure ratio. The expansion wave has been approximated by a negative shock. 7 = 1.4. 
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Figure 3>- Plot of average skin-friction coefficient against Reynolds number. 
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Figure 4.- Compressibility corrections as a function of shock pressure 
ratio for shock-tube flow. .7 = 1.4; cd = 0.8. 
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Figure 5.- Variation of compressible linear attenuation functions with 
shock pressure ratio. 7 = 1 . 4 j co = 0 . 8 . 
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(b) Region 3 . 


Figure 5*- Concluded 
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(a) Laminar flow; m = 1 . 

Figure 6 .- Plot of total compressible linear attenuation functions 
against shock pressure ratio. 7 =1.4; co = 0 . 8 . 
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Figure 6.- Continued. 
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(c) Turbulent flow; m = 13. 
Figure 6.- Concluded. 
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Figure 7«- Ratio of linear attenuation effects generated along forvard 
running characteristic to total attenuation effects of region 3 as 
a function of shock pressure ratio. 
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Distance from diaphragm , x 

Figure 8.- General wave diagram for nonlinear analysis in region |3. 
The short-dashed line is an extrapolation of the shock-wave path 
and the dotted lines are extrapolations of the characteristics. 
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(a) Details for computations of region I. 



(t>) Details for transitional flow. 

Figure 9 .- Detailed wave diagrams for nonlinear analysis in region 0. 
Short-dashed lines indicate extrapolations. 
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Figure 10.- Regional transition Reynolds number per foot of 
travel per atmosphere of ambient pressure plotted agains 
sure ratio. Circled symbols denote regions. 
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(a) Turbulent flow in a; laminar flow in (3. 


Figure 11.- Ratio of nonlinear to linear attenuation plotted against shock-wave position for 

^2. = 4.0 and = o.Ol x 10 6 . 
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(b) Turbulent flow in both a and (3. 


Figure 11.- Concluded. 
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Figure 12.- Ratio of nonlinear attenuation to linear attenuation for 
region a turbulent and. region p laminar plotted against shock pres- 
sure ratio. The shock wave is located at a distance Z/D from the 
diaphragm station. 
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Figure 12.- Concluded 
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Figure 13 •- Ratio of nonlinear attenuation to linear attenuation for 
both regions a and (3 turbulent plotted against shock pressure ratio. 
The shock wave is located at a distance l / D from the diaphragm 
station. 


Ratio of nonlinear to linear attenuation, 


NACA TN 4347 


77 





Figure 13*- Concluded 
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Figure l4. - Ratio of nonlinear attenuation to linear attenuation for 
region a turbulent and. region p transitional plotted against shock 
pressure ratio. R* = 1.25 x 106 . The shock wave is located at a 
distance Z/D from the diaphragm station. 
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Figure l4.~ Concluded 


Ratio of nonlinear to linear attenuation, • 


80 


NACA TN 4547 







0 2.4 6 8 10 12 14 16 18 20 

Shock pressure ratio, P^o/Pqq 


Figure 15 •- Ratio of nonlinear attenuation to linear attenuation for 
region a turbulent and region 0 transitional plotted against shock 
pressure ratio. R* = 2.50 X 10^. The shock wave is located at a 
distance l / D from the diaphragm station. 
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Figure 15 • - Concluded 
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Figure l6.- Schematic diagram of experimental setup for obtaining attenuation data. Distances 

are in feet. 



Figure 17.- Photograph of experimental setup shoving shock- tube low-pressure 

velocity- measuring stations. 
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Figure l8.- Photograph shoving details of velocity-measuring optical system. 
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Figure 19.- Schematic diagram of photomultiplier and output stages for shock velocity-measuring 

system. 
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Figure 20.- Plot of experimental averaged local shock pressure ratio plotted against shock 

distance from diaphragm. 
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Figure 21.- Theoretically predicted local shock pressure 

from diaphragm. 
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(b) ^ = 10.0. 


Figure 21.- Concluded 
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Figure 22.- Nonlinear theoretically predicted shock pressure ratios with transition plotted 

against shock distance from diaphragm. 
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Figure 22.- Concluded. 
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Figure 23.- Comparison of averaged experimental local shock pressure 

Pro 

ratios with various theoretical predictions. — — = 1.5. 
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Figure 24.- Comparison of averaged experimental local shock pressure 

ratios with various theoretical predictions. — ^ = I. 75 . Flagged 

-^00 

and unflagged symbols are used to distinguish between averages from 
different days. 
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Figure 25.- Comparison of averaged experimental local shock pressure 

Pgo 

ratios with various theoretical predictions, -t— =2.0. Flagged 

Poo 

and unflagged symbols are used to distinguish between averages from 
different days. 
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Figure 26.- Comparison of averaged experimental local shock pressure 


ratios with various theoretical predictions. 
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Figure 27.- Comparison of averaged experimental local shock pressure 
ratios with various theoretical predictions. —^2. =2.6. 
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Figure 28.- Comparison of averaged experimental local shock pressure 

J) Q - 

ratios with various theoretical predictions. = 3*0* Flagged 

-^00 

and unflagged symbols are used to distinguish between averages from 
different days. 
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Figure 29 .- Comparison of averaged experimental local shock pressure 

ratios with various theoretical predictions. = 3*6. Flagged 

-^00 

and unflagged symbols are used to distinguish between averages from 
different days. 
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Figure 30.- Comparison of averaged experimental local shock pressure 

ratios with various theoretical predictions. = b.O. The verti- 

-^00 

cal lines in (a) and (h) indicate the range of scatter in an averaged 
group when the scatter exceeds the symbol size. Flagged and un ag 
ged symbols are used to distinguish between averages from different 
days. 
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Figure 31 •- Comparison of averaged experimental local shock pressure 
ratios with various theoretical predictions. = 5.0. Flagged 

•^00 

and unflagged symbols are used to distinguish between averages from 
different days. 
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Figure 32.- Comparison of averaged experimental local shock pressure 

£62 = 6.0. 


ratios with various theoretical predictions. 


- 00 


102 


NACA TN 43A7 


l 7 , 1 

L?i 7 

“ nl 7j1 

NL 7 ,7 ° 

- — — fjL 7 »R — 2 .50> 
NL 7 ’j;R* =1 * 25> 

Experimental £ 

clO 6 

<10 6 

Lverag( 

53 




o 

o 

o 

. ,_a 



c 








-- 

— 






;a) T§ 

1 

01 X l! 

0 6 . 





' ~ . 





o 

0 











' ~ ~ ' 



— 

-ITT? 

i 


— 


(b) 

S?*o. 

02 x 1 

i 

0 6 . 

i 






























— 


(d) Ig 

- = 0.10 x 1( 

)*. 

1 



=. » — 


---- 

" " " 

























(e) T§” = 0.20 x ] 

i 1 ■ 


i 



















/ X* \ ^ 

(f) u a 

f-o . 

50 x 1 

0 6 . 



rr?ri 

-rri 

o o 

n O 


Q 2 4 6 8 10 12 M 16 IS 20 22 


Distance from diaphragm, ft 


Figure 33 


- Comparison of averaged experimental local shock pressure 
ratios with various theoretical predictions. = 7 - 0 - 
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Figure 34.- Comparison of averaged experimental local shock pressure 

Pq 

ratios with various theoretical predictions. — ^ =8.0. Flagged 

Pco 

and unflagged symbols are used to distinguish between averages from 
different days. 


104 


NACA TN 4347 





=4 

a 

% 

zrrl 







(b) 

?-• 

010 X 

10 6 . 






















'-■= 

— - v 

==== = 


= - - - 










- - - - 







W) iS 

f-o. 

050 x 

10 6 . 






- 




vigd 





















(•) ^=0.100 x 

10 6 . 






















SH 


















— 




200 x 

10 6 . 








2 4 6 8 10 12 14 16 18 20 

Distance from diaphragn, ft 


Figure 35 •- Comparison of averaged experimental local shock pressure 

Pr 

ratios with various theoretical predictions. — £2- = 9 . 0 . Flagged 

and unflagged symbols are used to distinguish between averages from 
different days. 
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Figure 36 .- Comparison of averaged experimental local shock pressure 

ratios with various theoretical predictions. — — = 10.0. The verti- 

-P 00 

cal lines in (a) and (b) show range of scatter in an averaged group 
when scatter exceeds symbol size. Flagged and unflagged symbols are 
used to distinguish between averages from different days. 
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